function h = circle2(x,y,r,color)
d = r*2;
px = x-r;
py = y-r;
rectangle('Position',[px py d d],'Curvature',[1,1],'FaceColor',color);
daspect([1,1,1])